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Neurons are excitable cells capable of sustaining high 
amplitude membrane potential oscillations called action 
potential (AP). There are two classes of repetitively firing 
neurons: Class 1 characterized by a continuously and 
arbitrarily low tunable firing frequency, or Class 2 char- 
acterized by a minimum firing frequency threshold [1]. 
The ionic channels and the morphology of the cells 
determine the characteristics of repetitive firing, such as 
the intrinsic period Pi, the firing threshold, or the dura- 
tion of the AP. For the study of phase locked modes and 
synchrony of neural populations a phenomenological 
understanding of neural activity suffices. In such cases, 
the phase response curve (PRC) rather than a detailed 
biophysical model of the neuron provides accurate and 
quick results. A PRC tabulates the transient changes in 
the firing period of a neuron due to an external perturba- 
tion, usually a brief current or synaptic conductance 
pulse [2,3]. The most significant effect of a perturbation 
occurs during the cycle that contains the perturbation 
and is quantified by the first order PRC, i.e. Fi((/)) = Pi/Pj 
- 1, where Pi is the transiently modified duration of the 



current cycle due to a perturbation applied at the stimu- 
lus time ts or phase (j) = ts/Pi (Figure lA). There are two 
types of PRCs: Type II PRC that has an almost sinusoidal 
shape and comparable sizes of the positive and negative 
lobes, and Type I PRC that has a disproportionate ratio 
of the two lobes (see Figure IB). When Type I neurons 
are part of neural network they can predominantly speed 
up (or only slow down) their rhythm, whereas Type II 
neurons synchronize more readily due to their ability to 
both slow down and speed up. While it is clear that every 
PRCs are bimodal, the purpose of this study is to investi- 
gate what parameters determine their shape. We found 
that the phase space bifurcation structure determines not 
only the PRC's type but also how easily the shape of the 
PRC can be tuned such that a "typical" Type I PRC can 
be smoothly converted into a Type II PRC. We found that 
nonlinear oscillators characterized both by saddle-node of 
invariant circle (SNIC) and Andronov-Hopf (HB) bifurca- 
tions, such as Morris-Lecar (ML) model neuron, the PRC 
can be smoothly adjusted from a pure Type I to a Type II. 
For example, by slowly increasing a bias current, ML 
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Figure 1 (A) Intrinsic oscillations (continuous line) and transiently perturbed period (dashed line). (B) A Type I PRC is bimodal with one large 
lobe of amplitude M and the other with a disproportionately smaller amplitude m. 



* Correspondence: oprisans@cofc.edu 

Department of Physics and Astronomy, College of Charleston, Charleston, SC 
29424, USA 

© 2014 Oprisan et al; licensee BiolVled Central ttd. This is an Open Access article distributed under the terms of the Creative Commons 
Attribution License (httpy/creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in 
any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http:// 
creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated. 



(3 Bion/led Central 



Oprisan et at. BMC Neuroscience 2014, 15(Suppl 1):P75 
http://www.biomedcentral.com/1471-2202/15/Sl/P75 



Page 2 of 2 



model reaches a SNIC bifurcation and produces a typical 
Type I PRC approximately captured by the generic infini- 
tesimal PRC iPRCi((/)) = 1 - cos{(f>). Further increasing the 
bias current brings the model closer to a HB that produces 
a typical Type II infinitesimal PRC iPRC2((/>) = sin{4> + a), 
with a a constant phase. For any intermediate value of the 
bias current the actual shape of the infinitesimal PRC is a 
linear combination of iPRCi and iPRC2. 
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